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ABSTRACT 

We report one of the most accurate measurements of the three-dimensional large-scale 
galaxy power spectrum achieved to date, using 56,159 redshifts of bright emission-line 
galaxies at effective redshift z « 0.6 from the WiggleZ Dark Energy Survey at the 
Anglo- Australian Telescope. We describe in detail how we construct the survey selec- 
tion function allowing for the varying target completeness and redshift completeness. 
We measure the total power with an accuracy of approximately 5% in wavenumber 
bands of Ak — 0.01 h Mpc -1 . A model power spectrum including non-linear correc- 
tions, combined with a linear galaxy bias factor and a simple model for redshift-space 
distortions, provides a good fit to our data for scales k < 0.4 h Mpc -1 . The large-scale 
shape of the power spectrum is consistent with the best-fitting matter and baryon 
densities determined by observations of the Cosmic Microwave Background radiation. 
^ ■ By splitting the power spectrum measurement as a function of tangential and radial 

wavenumbers we delineate the characteristic imprint of peculiar velocities. We use 
these to determine the growth rate of structure as a function of redshift in the range 
0.4 < z < 0.8, including a data point at z = 0.78 with an accuracy of 20%. Our growth 
rate measurements are a close match to the self-consistent prediction of the ACDM 
model. The WiggleZ Survey data will allow a wide range of investigations into the 
cosmological model, cosmic expansion and growth history, topology of cosmic struc- 
ture, and Gaussianity of the initial conditions. Our calculation of the survey selection 
function will be released at a future date via our website wigglez . swin . edu . au. 
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1 INTRODUCTION 

The pattern of density fluctuations in the low-redshift Uni- 
verse results from the physical processes which govern the 
evolution of matter perturbations after the Big Bang. In the 
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early Universe, the primordial spectrum of fluctuations cre- 
ated by inflation is processed before recombination in a man- 
ner depending on the physical matter density, baryon frac- 
tion and massive neutrino fraction (e.g. Bond & Efstathiou 
1984; Bardeen et al. 1986; Holtzman 1989; Hu & Sugiyama 
1996; Eisenstein & Hu 1998). After recombination, pertur- 
bations of all scales are amplified by gravity at an identical 



2 Blake et al. 



rate whilst linear theory applies. This growth rate depends 
on the matter and dark energy components which drive the 
cosmic expansion (e.g. Heath 1977; Hamilton 2001; Linder 
& Jenkins 2003; Percival 2005). The growth of fluctuations 
enters a non-linear regime at progressively larger scales at 
lower redshifts: in today's Universe, only perturbations with 
Fourier wavescales k < 0.1 h Mpc -1 evolve linearly to a good 
approximation (e.g. Smith et al. 2003; Jeong & Komatsu 
2006; McDonald 2007). 

The clustering pattern of galaxies at different redshifts 
is related to the underlying density fluctuations and may be 
used to test this model of structure formation. The shape of 
the clustering power spectrum - the relative amplitudes of 
large-scale and small-scale modes - depends on the composi- 
tion of the early Universe and may be used to extract infor- 
mation about the matter and baryon fractions (e.g. Tegmark 
et al. 2004a; Cole et al. 2005; Percival et al. 2007b). The 
amplitude of the clustering power spectrum as a function of 
redshift, together with the pattern of redshift-space distor- 
tions induced by galaxy peculiar velocities, can be used to 
measure the growth rate of structure (e.g. Hamilton 1992; 
Hawkins et al. 2003; Guzzo et al. 2008; Percival & White 
2009) . Higher-order or topological descriptors of the density 
field, such as the bispectrum or genus, can be applied to 
test whether the initial conditions are consistent with scale- 
invariant Gaussian random perturbations generated by in- 
flation (e.g. Gott, Dickinson & Melott 1986; Fry & Scherrer 
1994; Sefusatti & Komatsu 2007; James, Lewis & Colless 
2007). 

The interpretation of the shape and amplitude of the 
galaxy power spectrum is complicated by several factors. 
Firstly, the manner in which galaxies trace the density field - 
the "galaxy bias" - is in general a complex function of scale, 
dark matter halo mass, galaxy type and redshift (Dekel & 
Lahav 1999; Tegmark & Bromley 1999; Wild et al. 2005; 
Conway et al. 2005; Percival et al. 2007a; Smith, Scocci- 
marro & Sheth 2007; Cresswell & Percival 2009). However, 
the bias of galaxy fluctuations on sufficiently large scales 
(k < 0.1 h Mpc -1 at z = 0) appears to be well-described by 
a simple constant of proportionality whose value depends 
on galaxy type and luminosity, or more fundamentally dark 
matter halo mass (Peacock & Dodds 1994; Scherrer & Wein- 
berg 1998; Verde et al. 2002). Secondly, small-scale density 
perturbations eventually begin to evolve in a non-linear fash- 
ion requiring more complex modelling techniques such as 
higher-order perturbation theory or numerical iV-body sim- 
ulations (Smith et al. 2003; Jeong & Komatsu 2006; McDon- 
ald 2007). Thirdly, there is a practical challenge of acquiring 
galaxy survey data across a "fair sample" of the Universe 
(Tegmark 1997). For the large-scale linear modes of cluster- 
ing, which provide the most robust link to underlying theory, 
this sample must map a volume of the order 1 Gpc 3 using 
of the order 10 5 galaxies. These demands require multi-year 
campaigns with ground-based telescopes utilizing hundreds 
of clear nights (Glazebrook & Blake 2005). 

Despite these challenges, a series of galaxy redshift sur- 
veys have been undertaken to provide such datasets at red- 
shifts z < 0.5. The state-of-the-art projects which have 
mapped the "local" (z ~ 0.1) Universe are the 2-degree 
Field Galaxy Redshift Survey (2dFGRS; Colless et al. 2001) 
and the Sloan Digital Sky Survey (SDSS; York et al. 2000). 
The 2dFGRS obtained redshifts for 2 x 10 5 galaxies covering 



1500 deg 2 in the period between 1997 and 2002. The " main 
spectroscopic survey of the SDSS gathered 8 x 10 galaxy 
redshifts over 8000 deg 2 between the years 2000 and 2005. 
The SDSS project also included observations of 1 x 10 5 Lumi- 
nous Red Galaxies (LRGs) reaching up to a redshift z — 0.5 
(Eisenstein et al. 2001). 

These datasets have provided a rich source of informa- 
tion about the clustering of galaxies. For example, power 
spectra have been extracted for the 2dFGRS by Percival et 
al. (2001) and Cole et al. (2005); for the SDSS "main" galaxy 
sample by Pope et al. (2004), Tegmark et al. (2004a) and 
Percival et al. (2007a); and for the LRGs by Eisenstein et al. 
(2005), Huetsi (2006), Tegmark et al. (2006) and Percival et 
al. (2007b). Analysis of these surveys, in combination with 
the Cosmic Microwave Background fluctuations, has con- 
firmed that we inhabit a low-density Universe where matter 
today provides only 25 — 30% of the total energy governing 
the large-scale dynamics, with the rest located in a mys- 
terious "dark energy" component. In addition the baryonic 
fraction of the matter is only 15 — 20%, with the remainder 
composed of non-baryonic, cold particles whose nature is 
currently unknown (e.g. Percival et al. 2002; Tegmark et al. 
2004b; Tegmark et al. 2006; Komatsu et al. 2009). The clus- 
tering pattern is also sensitive to the presence of hot dark 
matter such as massive neutrinos, which comprise a small 
fraction of the energy budget (Elgaroy et al. 2002; Seljak et 
al. 2005). 

These galaxy surveys also describe how the underlying 
density fluctuations are modulated by galaxy bias (Verde et 
al. 2002; Wild et al. 2005; Conway et al. 2005; Percival et al. 
2007a; Cresswell & Percival 2009). In this context the com- 
parison of power spectrum measurements from the 2dFGRS 
and SDSS, which targeted galaxy populations selected in 
blue and red optical wavebands respectively, is of particular 
interest. When the differing galaxy types in these surveys 
are assigned linear bias factors, the resulting model fits to 
the linear-regime power spectra produce best-fitting mat- 
ter densities which are inconsistent at the statistical level 
of 2a. Careful treatment of scale-dependent and luminosity- 
dependent galaxy bias can potentially explain this discrep- 
ancy (Percival et al. 2007a; Sanchez & Cole 2008). 

There are strong motivations for extending these large- 
scale structure measurements to higher redshifts (z > 0.5). 
Firstly, the growth of structure implies that the linear regime 
of evolving perturbations extends to smaller scales at higher 
redshifts, enabling cleaner and more accurate model fits. 
Secondly, the shape of the survey cone allows access to signif- 
icantly greater cosmic volumes at higher redshift, enabling 
more accurate determinations of the large-scale power spec- 
trum amplitude. Thirdly, baryon oscillations in galaxy power 
spectra at different redshifts may be used as a standard ruler 
to extract the cosmic distance-redshift relation and infer the 
properties of dark energy (Blake & Glazebrook 2003; Seo & 
Eisenstein 2003; Hu & Haiman 2003). Fourthly, measure- 
ments of the growth of cosmic structure as a function of 
redshift increases our ability to discriminate between dark 
energy models including modifications to Einstein's theory 
of gravity (Guzzo et al. 2008; Wang 2008; White, Song & 
Percival 2009). 

Our current tools for probing the matter power spec- 
trum at redshifts z > 0.5 are limited. The clustering of 
high-redshift quasars has been studied by the 2dF Quasar 
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Survey (Outram et al. 2003) and the SDSS (Ross et al. 
2009) but the scarcity of QSOs implies that the large-scale 
clustering measurements are strongly limited by shot noise. 
Photometric redshifts from imaging surveys have been used 
to study the projected clustering pattern in redshift slices 
(Blake et al. 2007; Padmanabhan et al. 2007). However, this 
approach loses the information from small-scale radial clus- 
tering modes (Blake & Bridle 2005) and in particular pre- 
vents the extraction of the patterns of peculiar velocities, 
which are swamped by photometric redshift errors. Alterna- 
tively, fluctuations in the Lyman-a; forest absorption spec- 
trum on the sight lines to bright quasars have been used to 
infer the amplitude of small-scale clustering fluctuations in 
the high-redshift Universe (Croft et al. 2002; McDonald et 
al. 2005; McDonald et al. 2006). However, this method is 
potentially susceptible to systematic modelling errors and 
is only applicable at redshifts z > 2.3 where the Lyman-a 
absorption lines pass into optical wavebands. 

The WiggleZ Dark Energy Survey at the Anglo- 
Australian Telescope (Drinkwater et al. 2010) will provide 
the next step forwards in large-scale spectroscopic galaxy 
redshift surveys, mapping a cosmic volume of the order 1 
Gpc 3 over the redshift range z < 1. The survey, which be- 
gan in August 2006 and is scheduled to finish in July 2010, 
is obtaining of the order 200,000 redshifts for UV-selected 
emission-line galaxies covering of the order 1000 deg 2 of 
equatorial sky. The principal scientific goal is to measure 
baryon oscillations in the galaxy power spectrum in redshift 
bins up to z = 1 and provide a robust measurement of the 
dark energy model. The dataset will also trace the density 
field over unprecedented cosmic volumes at z > 0.5, pro- 
viding a sample comparable to the SDSS LRG catalogue at 
z < 0.5. Moreover, the spatial overlap between the WiggleZ 
and LRG catalogues in the redshift range 0.3 < z < 0.5 will 
allow careful studies of the systematic effects of galaxy bias 
on power spectrum estimation. 

This paper presents a determination of the current Wig- 
gleZ survey selection function and galaxy power spectrum, 
using a dataset comprising of the order 25% of the final 
survey observations. The selection function, which describes 
the angular and radial survey coverage in the absence of 
clustering, is complicated by the relatively high level of in- 
completeness in the survey affecting both the parent target 
catalogues and the spectroscopic follow-up observations (al- 
though this latter type of incompleteness will decrease as the 
survey progresses). However, we demonstrate that despite 
these complications the galaxy clustering power spectrum 
may be successfully extracted and already provides accu- 
rate tests of the cosmological model that rival lower-redshift 
surveys. 

The structure of this paper is as follows: in Section 2 
we present a detailed account of the survey selection func- 
tion including the coverage masks, completeness of our UV 
imaging catalogues, variations in redshift completeness, red- 
shift distribution as a function of sky position, and redshift 
blunder rate. In Section 3 we describe our power spectrum 
calculation and its correction for redshift blunders. We com- 
pare the predictions of cosmological models to the resulting 
power spectra in Section 4 and itemize our conclusions in 
Section 5. We note that we construct our selection function 
for a fiducial flat ACDM cosmological model with matter 
density O m = 0.3. 



2 WIGGLEZ SURVEY SELECTION 
FUNCTION 

An overview of the WiggleZ Survey observing strategy and 
galaxy selection criteria is presented by Blake et al. (2009) 
and Drinkwater et al. (2010) and summarized in Table 1. 
Briefly, targets are chosen by a joint selection in UV and 
optical wavebands, using observations by the Galaxy Evo- 
lution Explorer satellite (GALEX) matched with ground- 
based optical imaging. A series of magnitude and colour 
cuts is used to preferentially select high-redshift star-forming 
galaxies with bright emission lines, which are observed us- 
ing the AAOmega multi-object spectrograph at the Anglo- 
Australian Telescope (AAT). 

The survey selection function W(x) expresses the ex- 
pected mean density of galaxies with spectroscopic redshifts 
at position x, given the angular and luminosity survey se- 
lection criteria. An accurate determination of the selection 
function is essential in order to estimate the power spec- 
trum, which describes the amplitude of fluctuations relative 
to the mean density. 

Determination of the WiggleZ survey selection function 
is complicated by several factors, which are discussed in the 
following sub-sections: 

• The boundaries of each survey region are determined by 
the overlapping network of UV and optical imaging coverage 
[see Section 2.1]. 

• The magnitude cuts used to select galaxies from the in- 
put imaging lie close to the faint completeness thresholds of 
those surveys. For example, the completeness of the GALEX 
UV data at the survey flux threshold varies with the amount 
of Galactic foreground dust, inducing a variation in target 
density with angular position [see Section 2.2]. 

• The spectroscopic redshift completeness of each 
AAOmega pointing (i.e. the fraction of observed spectra pro- 
ducing successful redshifts) varies considerably with observ- 
ing conditions such as astronomical seeing and cloud cover, 
inducing a variation in the density of successful redshifts 
with angular position [see Section 2.3]. 

• The program of GALEX imaging has proceeded simul- 
taneously with the spectroscopic follow-up at the AAT. The 
expansion of the angular mask of the survey with time must 
be tracked in order to properly model the angular density 
of redshifts [see Section 2.3]. 

• The spectroscopic redshift completeness of each 
AAOmega pointing varies systematically across the spec- 
trograph field-of-view, decreasing towards the edges where 
acquisition errors are amplified [see Section 2.4]. 

• Whilst the survey is unfinished, the optical magnitude 
distribution of observed galaxies varies with position on the 
sky owing to the target prioritization scheme used for the 
observations (Drinkwater et al. 2010). This implies that the 
galaxy redshift distribution also varies with angular position 
[see Section 2.5]. 

• A fraction of the assigned galaxy redshifts are "blun- 
ders" resulting from emission-line mis-identifications. The 
rate of blunders depends on redshift, and hence distorts the 
true radial galaxy distribution [see Section 2.7]. 

2.1 Coverage masks 

This part of the selection function establishes a (0,1) binary 
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Table 1. Summary of the WiggleZ Survey galaxy selection criteria. 



UV magnitude cut NUV < 22.8 

Optical magnitude cut 20 < r < 22.5 

UV colour cut FUV dropout or FUV - NUV > 1 

Optical-UV colour cut -0.5 < NUV - r < 2 

Optical colour cut If g < 22.5 and i < 21.5, then (r — i) < (g — r — 0.1) and (r — i) < 0.4 



angular coverage mask indicating the availability of input 
targets. 

The boundaries of each WiggleZ survey region are de- 
fined by the coverage maps of the UV and optical imaging 
data within that region. The UV data consists of a series 
of pointings of the GALEX satellite. The GALEX field-of- 
view is a circle of radius 0.6 deg; we select sources from a 
slightly smaller radius of 0.55 deg due to concerns about the 
GALEX photometry at the edges of the field (Morrissey et 
al. 2007). For the survey regions analyzed in this paper, the 
optical coverage map corresponds to the 4th Data Release 
of the SDSS. 

The GALEX source catalogues contain small regions of 
bad data corresponding to scattered light from bright stars 
adjacent to the field-of-view. We inspect each GALEX field 
by eye for the presence of these artefacts and define rectan- 
gular masks to remove this bad data. These masks encom- 
pass a negligible fraction (< 1%) of the survey area. 

2.2 Variation of parent density with dust and 
exposure time 

This part of the selection function modulates the angular cov- 
erage map in accordance with incompleteness in the parent 
imaging catalogues. 

The faint UV magnitude threshold for WiggleZ survey 
selection, NUV = 22.8, is comparable to the magnitude 
completeness limit of the GALEX Medium Imaging Survey 
(MIS) which provides the targets. The incompleteness of the 
WiggleZ target catalogue is therefore significant and must 
be modelled in our selection function, because it determines 
the baseline target density on the sky in the absence of clus- 
tering. We calibrated this incompleteness factor as a func- 
tion of foreground Galactic dust (quantified by the value of 
Eb-v) and GALEX exposure time t cxp . 

The variations of Eb-v and i cxp across the survey 
regions analyzed in this paper are displayed in Figure 1. 
The distribution of exposure time is peaked in the range 
1400 < icxp < 1800 sec (the canonical MIS exposure time 
is 1500 sec, equivalent to one orbit of the GALEX satellite 
when overheads are included). However, our analysis also 
includes fields with lower exposure times 500 < t oxp < 1400 
sec for which we are still gathering data. In addition some 
fields have significantly larger values of t cxp if they have 
been observed as part of other GALEX projects. The distri- 
bution of dust extinction is broad, peaking at Eb-v ~ 0.04. 
Although we correct the UV magnitudes for dust extinction 
(by an amount Anuv ~ 8 Eb-v ~ 0.3, see Drinkwater et 
al. 2010), the proximity of our observations to the MIS mag- 
nitude threshold induces a variation of completeness with 
Eb-v- 

We used the GALEX NUV (dust-corrected) differen- 
tial source counts (i.e., counts of galaxies in NUV magni- 
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Figure 1. Histograms of the values of GALEX exposure time 
and Galactic dust extinction across the WiggleZ survey regions 
analyzed in this paper. 



tude bins) to calibrate the dependence of the WiggleZ survey 
completeness on Eb-v and i exp . For this analysis we created 
a GALEX-SDSS matched galaxy sample for which we did 
not impose any of the WiggleZ survey magnitude or colour 
cuts. This is because the WiggleZ selection cuts reduce the 
target density by a factor of approximately 6, greatly in- 
creasing the noise in these measurements. We matched the 
GALEX NUV catalogues with the SDSS data (using a tol- 
erance of 2.5 arcsec) and removed objects flagged as stars in 
SDSS. The presence of these stars would distort the galaxy 
number counts and induce spurious correlations with Eb-v 
owing to the non-uniform stellar distribution. 

We constructed the source count by dividing the survey 
coverage map into small pixels of size 0.1 x 0.1 deg and 
assigning each pixel a mean value of E B -v and i cxp . We 
then added up the pixel source counts in bins of Eb-v and 

^cx P • 

Figure 2 illustrates the source count in bins of exposure 
time, restricting the analysis to pixels with low dust extinc- 
tion (Eb-v < 0.04) in order to isolate the variation with 
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counts as a function of de-reddened magnitude m = NUV 
is a power-law modulated by an incompleteness function: 



NUV (de-reddened) 

Figure 2. Differential source counts of GALEX-SDSS galaxy 
matches in bins of GALEX exposure time. We restrict the analysis 
to regions of the survey with low dust extinction Eb—v < 0.04. 
We note that the uppermost curve with the caption "6000 < 
cxp < 99999" corresponds to the "fiducial" source count in the 
limit of low dust extinction and high GALEX exposure time, as 
discussed in the text. 
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Figure 3. Differential source counts of GALEX-SDSS galaxy 
matches in bins of dust extinction. We restrict the analysis to 
regions of the survey with GALEX exposure times 1400 < f. eX p < 
1800 sec. 



exposure time. We expect to see the number counts rising 
with increasing magnitude (in the classic Euclidean regime, 
d(log 10 N)/dm = 0.6) and tailing off at faint magnitudes 
when incompleteness becomes significant. Indeed, Figure 2 
illustrates that the counts in different exposure time bins 
agree well at bright magnitudes, and the completeness limit 
grows fainter with increasing exposure time. More detailed 
fits are presented below. 

Figure 3 displays the counts in bins of Eb-v for 
GALEX observations in the exposure time range 1400 < 
£ CX p < 1800 sec, in order to isolate the variation with dust 
extinction. Again, the agreement is good at bright magni- 
tudes (i.e. there is no significant residual dependence of the 
number counts on dust extinction) and the incompleteness 
at faint magnitudes increases with the value of Eb-v- 

We found that a good fitting formula for these source 



dN dN 

— — (m) = — — (m) x C(u,a,m) 
dm dm 



where 
dN 

C(fj,, a, m) 
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where erf (a;) is the error function: 



erf (a;) = -^= / exp (— t 2 ) dt 



(1) 

(2) 
(3) 

(4) 



We first fitted Equation 1 to the source counts for 
survey regions in the ranges Eb-v < 0.04 and t cxp > 
6000 sec, which we designated the "fiducial" source count 
Co ([to , o"o , m) in the limit of low dust extinction and high 
GALEX exposure time. The apparent incompleteness in this 
measurement is not due to GALEX observations but to 
the limiting magnitude threshold of the SDSS data with 
which we match the GALEX sources. We fitted the model 
to the magnitude range 20.0 < NUV < 22.8 (motivated 
by the faint flux threshold for WiggleZ target selection). 
The best-fitting parameters for the fiducial source count are 
P = 0.625±0.013, (i = 22.98±0.05 and a = 1.1±0.1, which 
we assumed for the rest of this analysis. The errors in each 
parameter are quoted after marginalizing over the remaining 
parameters. Our derived number counts are consistent with 
Xu et al. (2005), who show that a model incorporating lu- 
minosity evolution provides a good fit to this data. We note 
that the overall amplitude of the number counts, parame- 
terized by a in Equation 2, is not required for determining 
the incompleteness, as explained below. 

We then defined the incompleteness in the GALEX cat- 
alogues as a function of Eb-v and t cxp , relative to these 
fiducial counts, by fitting the model: 

dN dN . . 

— — (m) = — — (m) x Go (/io, <jo, m) x C(u,a,m) 
dm dm 



(5) 



By construction, C = 1 in the limit of low dust extinction 
and high exposure time. We found that the best-fitting value 
of a does not vary significantly with dust or exposure time. 
Therefore we fixed a — 0.578 such that the completeness 
function depended on just one parameter, fi. We then fitted 
the value of /i as a function of both Eb-v and t axp ; the 
results are displayed in Figure 4, together with the error in [i 
obtained by marginalizing over the other fitted parameters. 

We can compare these results to theoretical expecta- 
tions. The variation of the completeness limit /x with £ e x P 
(at fixed Eb-v) agrees well with the theoretical expecta- 
tion for background-limited imaging of fixed signal-to-noise, 
in which a doubling of exposure time equates to approxi- 
mately 0.4 magnitudes of survey depth. The variation of /i 
with Eb-v (at fixed t C x P ) should follow the dust-extinction 
law A/i w 8AEb-v- The observed slope in Figure 4 is in 
fact a little steeper, owing to the influence of the additional 
factor Co(m) in Equation 5. 

The completeness of the GALEX-SDSS parent cata- 
logue for the UV flux threshold mo = 22.8 can be deter- 
mined for a given p(t eX p, Eb-v) by evaluating 
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Figure 4. Dependence of the completeness parameter fi on dust 
extinction, in bins of exposure time. The measurements of u 
are of highest quality for the dominant exposure time range 
1400 < t CX p < 1800 sec. The error range in p, is asymmetric about 
the best-fitting value because the model becomes increasingly in- 
sensitive to /i as the value of /i increases. The measurement for 
each dust bin is plotted at the median value of Eb-v for that 
bin. 
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Figure 6. Comparison of the observed and modelled target den- 
sity for the WiggleZ catalogue in bins of dust extinction for the 
GALEX exposure time range 1400 < t exp < 1800 sec. The av- 
erage value of Eg_v for each dust bin is displayed next to the 
corresponding data point. Poisson error bars, which are likely to 
be a mild under-estimate of the true error, are shown for the 
measured density. 




20.5 21 21.5 

NUV (de-reddened) 

Figure 5. Fraction of GALEX-SDSS matches that are selected 
as WiggleZ targets as a function of GALEX NUV magnitude. 
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In order to model the WiggleZ survey target density 
from these completeness measurements we must now allow 
for the fraction of GALEX-SDSS matched sources that are 
selected as WiggleZ targets as a function of NUV magni- 
tude. This function is plotted in Figure 5. WiggleZ targets 
are preferentially faint in NUV owing to the colour selec- 
tion cut -0.5 < NUV - r < 2. We weighted Equation 6 
with this function in order to determine the WiggleZ sur- 
vey angular completeness map. We defined this map using 
a uniform grid of pixels in right ascension and declination. 

In Figure 6 we compare the WiggleZ target densities as 
a function of dust extinction predicted by our number-counts 
modelling with the densities observed in the catalogue. We 
plot the percentage difference of the measured density from 



the predicted density. The good match indicates that our 
selection function model is successful. An example parent 
catalogue completeness map for the 9-hr survey region is 
displayed as a panel in Figure 7. 

We note that the variation in the angular density of tar- 
gets is dominated by incompleteness in the UV imaging data 
rather than in the optical imaging data. Although the faint 
optical magnitude threshold for WiggleZ selection (r = 22.5) 
lies at the completeness limit of the SDSS, the NUV — r 
colour selection cut implies that the median r-band mag- 
nitude of WiggleZ targets is r ~ 21.5 and the variation of 
SDSS completeness with Eb-v is negligible. We also stud- 
ied the variation of completeness with the local astronomical 
seeing in the SDSS images, and found no significant effect 
on the target densities. 



2.3 Redshift completeness map 

This part of the selection function samples the density map 
of the parent target catalogue with the pattern of spectro- 
scopic follow-up observations. 

The spectroscopic observations of the WiggleZ survey 
are defined by a series of AAT pointing centres across each 
survey region. The pointing centres for each region are de- 
termined prior to each survey observing run based upon the 
available distribution of targets at the time, using a "simu- 
lated annealing" algorithm. Each patch of sky must be ob- 
served three to four times on average to build up the WiggleZ 
redshift catalogue, therefore we obtain a series of overlap- 
ping pointings. Each individual pointing results in a fraction 
of successful redshifts which varies considerably with astro- 
nomical seeing, airmass and cloud cover from under 40% in 
poor conditions (> 2.5 arcsec seeing or airmass sec(z) > 1.5 
or > 2/8ths cloud) to over 80% in good conditions (< 1.5 
arcsec seeing and sec(z) < 1.4 and no cloud). Furthermore, 
the fraction of successful redshifts in any pointing varies 
across the field plate of the 2dF spectrograph, as discussed in 
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Figure 7. Panels showing for the WiggleZ 9-hr region: (i) the angular completeness map of the parent catalogue, (ii) a realization of the 
parent catalogue sampled from this map, (iii) the sequence of 2dF field centres, (iv) a realization of the redshift catalogue. Note in panel 
(iii) that in some pointings of the survey a field radius of 0.7 deg rather than 1.0 deg was used, due to concerns over radial-dependent 
redshift completeness within 2dF pointings. 



more detail below. Therefore, as the survey is still partially 
complete, the redshift completeness map is a complicated 
function of position on the sky. 

We determined the redshift completeness function in 
each survey region using Monte Carlo realizations of the ob- 
servations. The steps in the construction of a realization are 
illustrated in Figure 7. For each region we started with the 
density map as a function of dust and exposure time created 
by the process described in Section 2.2 [Figure 7, panel (i)]. 
We created a Monte Carlo realization of this density map, 
containing the same number of galaxies as the real target 
catalogue, by first generating a uniform random catalogue 
and then excluding points based on a probability equal to 
the local incompleteness [Figure 7, panel (ii)]. The spectro- 
scopic observations consist of a series of overlapping AAT 
pointings [Figure 7, panel (iii)]. Each pointing was applied 
to our realization by laying down a field circle at the cor- 
rect angular position and randomly assigning the appropri- 
ate number of observed sources to spectroscopic fibres. 

We note that the available target area has grown over 
the duration of the survey because the GALEX imaging ob- 
servations are proceeding simultaneously with the redshift 
follow-up. In the Monte Carlo simulations it is therefore nec- 
essary to track the distribution of GALEX tiles available 
prior to each WiggleZ observing run and only assign tar- 
get sources to each random realization in areas that were 
available at the time of each pointing. 

For each telescope pointing in the random realization, 



a fraction of the observed sources are flagged with success- 
ful redshifts. The remainder are flagged with unsuccessful 
redshifts. The total number of successful redshifts for each 
pointing in the Monte Carlo realization is equal to that 
achieved in the corresponding survey observation. The prob- 
ability of a given random source being assigned a success- 
ful redshift varies across the field-of-view of the pointing in 
the manner described below in Section 2.4. The resulting 
random realization, constructed by applying the full set of 
pointings, is displayed in Figure 7, panel (iv). 

The WiggleZ survey observing strategy involves obtain- 
ing repeat spectra of unsuccessful redshifts observed in poor 
conditions. In addition, a small number of galaxies with ex- 
isting successful redshifts are re-observed in order to quan- 
tify the redshift blunder rate. Each of these categories of ob- 
servation are included when constructing the Monte Carlo 
realizations. 



2.4 Variation of redshift completeness across the 
spectrograph field-of-view 

The probability of obtaining a successful redshift in a Wig- 
gleZ survey pointing depends on the source position within 
the spectrograph field-of-view. We find that galaxies ob- 
served at the edges of the field plate have a reduced redshift 
completeness because these observations are more seriously 
affected by rotational mis-alignments in target acquisition, 
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Figure 8. Radial dependence of the rcdshift completeness in the 2-degree field for the six WiggleZ survey observing runs between August 
2007 and May 2008. In August 2007 the field of observation was restricted to radii < 0.7 degrees. We note that the completeness is 
normalized relative to a level of 1.0 at the centre of each field. 
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Figure 9. Dependence of the redshift completeness on position in the 2-degree field for the six WiggleZ survey observing runs between 
August 2007 and May 2008. In August 2007 the field of observation was restricted to radii < 0.7 degrees. We note that the completeness 
is normalized relative to a mean level of 1.0 averaged across each field. 
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which is performed using a limited number of guide-fibre 
bundles. 

Furthermore, we find that the pattern of redshift com- 
pleteness across the field-of-view varies between survey ob- 
serving runs as the spectrograph set-up is re-calibrated, but 
remains stable throughout an observing run. Although indi- 
vidual fibres can be placed in the field-of-view by the 2dF 
robot positioner with an accuracy of 20/xm (0.3 arcsec), sys- 
tematic errors can develop in the mapping of (x,y) posi- 
tion on the field plate to (R.A., Dec.) position on the sky. 
These errors are caused by minor discrepancies in the at- 
mospheric correction applied to target apparent positions, 
and (at times during the survey) by an instability in the 
modelling of opto-mechanical distortions between the prime 
focus corrector and the 2dF field plates. These additional 
errors imply that the redshift completeness variation is not 
radially symmetric. 

These effects are illustrated by Figures 8 and 9. In Fig- 
ure 8 we display the radial variation in redshift completeness 
for WiggleZ survey observing runs between August 2007 and 
May 2008. Figure 9 displays the 2D variation of redshift com- 
pleteness across the field plate for each of these observing 
runs, and we use these 2D maps to generate the probabili- 
ties of successful redshifts in the Monte Carlo realizations de- 
scribed above. In some observing runs, such as August 2007, 
we restricted our spectroscopic observations to the central 
0.7 deg radius of the 1-deg radius field-of-view of the 2dF 
spectrographs, because the drop in redshift completeness at 
the edges of the field was particularly significant. 
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Figure 10. The distribution of WiggleZ redshifts z in the five 
magnitude bands in which the targets are prioritized. Data is 
shown for the combined 9-hr, 11-hr and 15-hr survey regions for 
the redshift range 0.3 < z < 0.9 analyzed in this paper. Poisson 
error bars are shown for the points and the data is fitted with a 
sum of two Gaussian functions. 



for the five magnitude bands for the combined WiggleZ re- 
gions. We fit the redshift distributions with a sum of two 
Gaussian functions and sample the random redshifts from 
these smooth distributions. Future work will establish N(z) 
from measurements of the galaxy luminosity function in- 
stead. 



2.5 Radial selection function versus angular 
position 

This part of the selection function creates the appropriate 
radial distribution of galaxies depending on the distribution 
of r -band magnitudes in each 2dF pointing. 

The final step in the creation of the Monte Carlo real- 
izations is to assign a random redshift to each galaxy, thus 
establishing the selection function in the radial direction. 
Targets in the WiggleZ survey are prioritized for observa- 
tion in accordance with their SDSS r-band magnitudes such 
that the faintest galaxies are observed earliest in the se- 
quence (Drinkwater et al. 2010). In detail we use five prior- 
ity bands which divide the magnitude range 20.0 < r < 22.5 
into equal pieces. Given that there is a correlation between 
redshift and magnitude, this implies that the survey redshift 
distribution varies across the sky in a manner dependent on 
the density of redshifts obtained in a given area. 

We track this in our random catalogues by recording for 
each telescope pointing the number of successful redshifts 
obtained in each magnitude priority band, and assigning 
these magnitude identifications to sources in the random 
catalogues. We note that the weak dependence of redshift 
completeness on r-band magnitude is also absorbed into this 
step of the process. 

For each magnitude band we measured the redshift dis- 
tribution N(z) of successful redshifts using the existing Wig- 
gleZ spectroscopic data. In order to reduce the fluctuations 
due to cosmic variance, we combined the data for the 9- 
hr, 11-hr and 15-hr survey regions in this analysis. We used 
these probability distributions to assign redshifts to random 
sources on the basis of their magnitude. Figure 10 plots N(z) 



2.6 Final construction of the survey selection 
function 

The method described in the preceding sub-sections allows 
us to construct Monte Carlo realizations for each WiggleZ 
region incorporating the angular and radial variations of the 
selection function. We can accurately determine the full se- 
lection function function W(x) by stacking together many 
such realizations: we typically generate 10,000 realizations 
for each region. We pixelize the selection function in grid 
cells corresponding to our power spectrum measurement, as 
discussed in Section 3.1. When converting redshifts to dis- 
tances we use a fiducial flat ACDM cosmological model with 
matter density f2 m = 0.3. 

2.7 Redshift blunder rate 

The redshifts of galaxies in the WiggleZ survey are typi- 
cally based on identifications of emission lines; the signal- 
to-noise of the spectra is usually too low to permit de- 
tection of the galaxy continuum. The principal line used 
for redshift identification is the [Oil] doublet at rest-frame 
wavelength 3727A. This emission line lies in our observed 
spectral window 4700 — 9500A for the galaxy redshift range 
0.26 < 2 < 1.55. The redshift identification is confirmed 
for most galaxies by the additional presence of emission 
lines such as H/3 486lA, [OIII] 4959A, [OIII] 5007A, and 
Ha 6563A (Drinkwater et al. 2010). 

However, not all redshift identifications are based on 
multiple emission lines. Features redward of [Oil] progres- 
sively leave the observed spectral range with increasing red- 
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shift. The H/3 and Ha lines are observable for the ranges 
z < 0.95 and z < 0.45, respectively. At relatively high red- 
shifts the galaxy emission lines must be identified against a 
background of noisy sky emission lines. Despite these diffi- 
culties, we can gain some confidence in single-line redshifts 
based on [Oil] either through detection of the doublet, which 
is marginally possible with our spectral resolution for galax- 
ies lying at z > 0.8, or by eliminating other solutions by 
failure to detect [Oil] at lower wavelengths in cleaner parts 
of the spectrum. 

We assign quality flags from Q = 1 (lowest) to Q = 5 
(highest) for each WiggleZ redshift based on the confidence 
of our measurement. Redshifts with quality Q > 3 are con- 
sidered "reliable" and used in our analysis. Redshifts with 
quality Q > 4 are based on multiple emission lines and are 
very secure. Galaxies with redshifts based on noisy data or 
single emission lines are assigned Q = 3. The fraction of 
reliable redshifts with Q — 3 is approximately one-third. 

Some fraction of WiggleZ redshifts will be blunders. We 
distinguish two types of redshift blunder for the purposes of 
our analysis. Firstly, a galaxy emission line may be mis- 
identified as another, incorrect, emission line. In our power 
spectrum measurement this represents (approximately) a 
convolution of the galaxy density field whereby structures at 
a given redshift are coherently copied to a second redshift. 
Secondly, a night-sky emission line may be mis-identified 
as a galaxy emission line. As there are a large number of 
night-sky emission lines available for mis-identification, this 
effectively corresponds to a randomizing of the galaxy den- 
sity field through the addition of objects whose positions are 
uncorrelated with the underlying density. 

We studied the redshift blunder rate through a pro- 
gramme of repeat observations. In each survey pointing we 
assigned a small number of spectrograph fibres (typically 
3-5) to galaxies which have already been assigned redshifts 
with quality Q > 3. We define two redshifts as inconsistent 
if they differ by Az > 0.002 (the typical redshift error for 
our spectra is Az = 0.0005 or 100 km s _1 ). 

We find that pairs of repeat galaxy redshifts which both 
possess quality Q > 4 disagree in 2% of cases. Assuming that 
one of the pair of inconsistent values is the correct redshift, 
this implies that the blunder rate for the set of Q > 4 red- 
shifts is 1%. Pairs of repeat redshifts which both possess 
Q — 3 disagree in 31 ± 2% of cases. However, we can obtain 
a larger statistical sample for analysis if we consider pairs 
composed of Q — 3 and Q > 4 redshifts, supposing that 
the higher quality redshift is the correct value. Under this 
method we find that the blunder rate of Q = 3 redshifts 
is 17 ± 1%, in good agreement with the internal blunder 
rate for Q = 3 pairs. Given that approximately one-third 
of reliable redshifts are assigned Q = 3, the overall blunder 
rate for the WiggleZ survey is about 5%. However, we must 
carefully quantify the redshift blunders in more detail in or- 
der to obtain an unbiased measurement of the galaxy power 
spectrum. 

In Figure 11 we illustrate the fraction of redshift blun- 
ders resulting from the mis-identification of an emission line 
with a second, incorrect, emission line by plotting the distri- 
bution of values of (1 + Z\)/(l + £2) for inconsistent repeat 
redshifts composed of Q — 3 and Q > 4 pairs. The val- 
ues of {z\,Z2) are the redshifts of the Q — 3 and Q > 4 
spectra, respectively. This histogram reveals three signifi- 
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Figure 11. Distribution of values of (1 + z\)/(l + 22) for incon- 
sistent repeat redshifts derived from pairs of spectra with quality 
flags Q = 3 (redshift 21) and Q > 4 (redshift 22). The vertical 
lines indicate the ratios expected in the cases where H/3, [OIII] 
and Hc» are mis-identified as [Oil]. 



cant spikes corresponding to the mis-identification of H/3, 
[OIII] 5007A and Ha as [Oil]. Approximately 30% of red- 
shift blunders correspond to this type of mis-identification; 
the correct redshift in such cases is typically lower than the 
blunder redshift. 

Figure 12 plots the distribution of redshift blunders not 
contained in the three spikes in Figure 11. This type of blun- 
der, comprising about 70% of all blunders, corresponds to 
mis-identification of sky emission lines as [Oil]. In Figure 
12 we have also fitted a model for the redshift distribution 
N(z) of the form: 



N(z) <x 



exp 



(7) 



The best-fitting parameters are zo = 1.11, a = 1.24, 
ft = 6.46. We note that the distribution of blunders peaks 
at a significantly higher redshift than the distribution of 
correct redshifts shown in Figure 10. This implies that the 
blunder fraction varies significantly with redshift - this be- 
haviour is plotted in Blake et al. 2009, figure 6. In Section 
3.2 we model the effect of these types of redshift blunders 
on measurements of the galaxy power spectrum. 



3 POWER SPECTRUM ANALYSIS 

3.1 Power spectrum estimation methodology 

In this Section we summarize our method of power spectrum 
estimation, prior to presenting our analysis of the WiggleZ 
survey data in Section 3.3. Our power spectrum estimation 
is based on the optimal weighting scheme of Feldman, Kaiser 
& Peacock (1994) [FKP] (also see the discussions in Tadros 
& Efstathiou 1996; Hoyle et al. 2002). When converting red- 
shifts to distances we use a fiducial flat ACDM cosmological 
model with matter density f2 m = 0.3. We first enclosed the 
survey cone for the particular region and redshift interval 
within a cuboid of sides (L x , L y ,L z ) and gridded the galaxy 
catalogue in cells numbering (n x ,n y ,n z ) using nearest grid 
point (NGP) assignment to produce a distribution n{x) . The 
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Figure 12. The distribution of redshift blunders for the pairs in 
Figure 11 which do not lie in the three prominent spikes. The 
solid line is the best fit of the model described by Equation 7. 



cell dimensions were chosen such that the Nyquist frequen- 
cies in each direction (e.g. fcNyq.x = ^n x /L x ) exceeded the 
maximum frequency of measured p O we r Km ax by at least a 
factor of two (although we corrected for the effect of NGP 
assignment on the power spectrum as explained below). 

We then applied a Fast Fourier transform to the grid 
of data, weighting each pixel by a factor w(x) following the 
method of FKP: 



n(k) — S ' n(x)w(x) exp (ik.x) 

X 

where the weighting factor is given by 
1 



1 + W(x)N c n P 



(8) 



(9) 



n x n y n z is the 



(FKP, equation 2.3.4). In Equation 9, N< 
total number of grid cells, no = N/V is the mean number 
density of galaxies in units of h? Mpc" 3 , where N is the 
total number of galaxies and V = L x L y L z is the cuboid 
volume, and Po is a characteristic value of the power spec- 
trum at the Fourier wavescale of interest. The purpose of 
the weighting scheme of Equation 9 is give equal weight per 
volume where we are limited by cosmic variance, and equal 
weight per galaxy where we are limited by shot noise, re- 
sulting in an optimal measurement of the power spectrum 
amplitude P(fc) in the case where P(k) ~ Pq. 

In Equation 9, W(x) is proportional to the survey selec- 
tion function determined in Section 2, which describes the 
number of galaxies expected in each cell x in the absence of 
clustering. We assume here the normalization 



X) = 1. 



(10) 



We note that summations over cells x may be related to the 
equivalent integration over volume elements d 3 x by 



i-^W(f) = I J W (x)d 3 x. 



(11) 



We also apply a Fast Fourier transform to the survey selec- 
tion function: 



W(k) = > W(x)w(x) exp (ik.x) 



(12) 



The estimate of the galaxy power spectrum for wavescale k 
is then 

g \h(k) - NW(k)\ 2 ~ NJ2 s W(x)w 2 (x) 
cstW ~ N*N c J2 s W*(x)wHx) (6) 

This estimate of the power spectrum is biased by the process 
of NGP assignment ( Jing 2005) . In order to remove this bias 
we calculated the correction factor for each Fourier mode (by 
which the power spectrum estimate should be divided): 



Correction factor(fc) = 



J2 A H 2 (k>)P(k') 
P(k) 



(14) 



where m = (m x ,m y ,m z ) is a vector of integers, k 1 = 

{h xi ky , k z ) = {k X ~\~ fflx&Nyq.X J ky -\- TTLyk^y^y , k Z 

and 



H(k) 



( 7Tk x \ ( TTky \ / 7Tkz \ 



(15) 



In Equation 14, P(k) is the underlying model power spec- 
trum. Given that this is initially unknown, we proceed by an 
iterative approach: we assume a fiducial cosmological model, 
compute the correction factor, fit cosmological parameters 
to the power spectrum, re-calculate the correction factor, 
and then repeat the parameter fit. The magnitude of the 
correction is typically 2% at scale k ~ 0.2 h Mpc -1 . 

The spatially- varying selection function W(x) has two 
effects on the process of power spectrum estimation. Firstly 
the expectation value of Equation 13 is the underlying power 
spectrum P(k) convolved with the survey selection function: 



< Pest(fc) >= 



£ fc -,p(fcOivnfc-fc , ')i i 

N c Y /S W 2 (x)w 2 {x) 



(16) 



The numerator of Equation 16 is summed over the grid 
points k' in Fourier space for which the Fast Fourier Trans- 
form of W(x) is calculated, i.e. spaced by (Ak x , Ak y , Ak z ) = 
(2n/L x , 2-k/L v , 2n/Lz). For reasons of computing speed 
when fitting models, we re-cast this equation as a matrix 
multiplication in Fourier bins of width Afc = O.Ol/i Mpc -1 : 



< Post 



(17) 



We determine the convolution matrix My by evaluating the 
full sum of equation 16 for a set of unit vectors, e.g. for bin 



P(k) = 1 (hi .mm < |fc| < h 
= otherwise 



(18) 
(19) 



Secondly, the estimates of the power in different Fourier 
modes k become correlated. If we average the estimates of 
Equation 13 into bins in Fourier space, labelling the bins by 
i, the covariance between bins i and j is given by 

2 



< 5Pi SPi >= j^rp J2 \ P ®$ - + S(k - k') 



NgN,, — 

fc.fe' 



(20) 



(FKP, equation 2.5.2) where k and fc' are pairs of Fourier 
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modes lying in bins i and j, P is the characteristic power 
spectrum amplitude in bins i and j defined below, and the 
functions Q(k) and S(k) are given by FKP equations 2.2.3 
and 2.2.5: 



Q{k) 
S(k) 



W 2 {x)w 2 (x) exp (ik.x) 
1 \^2 s W(x)w 2 {x)exp(ik.x) 



n N c 



Y, B W*{x)w\x) 



(21) 
(22) 



In deriving Equation 20 it is assumed that the power spec- 
trum factor P which appears is effectively constant over 
Fourier separations k — k' which produce correlated esti- 
mates. For our datasets the Fourier transform of the selec- 
tion function, W(k), is sufficiently compact around k = 
that this is a valid approximation. We evaluated Equation 
20 for each survey region by a direct summation over Fourier 
modes in the FFT grid. Equation 20 depends on the underly- 
ing power spectrum, which is initially unknown, in the same 
manner as equation 14. We again used an iterative approach 
whereby we initially used a default model power spectrum to 
make this calculation, and then replaced it using the fitted 
parameters. 

In order to facilitate comparison with other studies it is 
useful to take the limit of these equations in the case where 
the selection function is constant, i.e. W(x) = 1/N C . The 
power spectrum estimator of Equation 13 becomes 



Pest(fe) 



\n(k) — NW(k)\ 2 — 7V 
TV 2 



(23) 



and the covariance matrix in Equation 20 reduces to a di- 
agonal matrix with entries: 



iVg V no/ 



(24) 



where is the number of Fourier modes lying in bin i. 
Equation 24 clarifies that there are two sources of error in 
an estimate of the power spectrum: cosmic variance and shot 
noise, represented by the two terms inside the bracket. 



3.2 Redshift blunder correction 

In this Section we calculate the distortion in the galaxy 
power spectrum created by the types of redshift blunder 
described in Section 2.7. In order to gain intuition we be- 
gin with a simple model using the "flat-sky approximation" 
which supposes that galaxies are scattered in position along 
a single axis of the cuboid (which we take as the x-axis). 
Defining S(x) as the galaxy overdensity in the cell at posi- 
tion x — (x,y,z), the galaxy number distribution is given 
by 



n(x) = NW(x)[l + S(x)] 



(25) 



where W(x) is the selection function normalized as above. 
We now suppose that a fraction / of galaxies are scattered 
in position along the a;-axis such that their final ^-position 
is drawn from a probability distribution V(x) (i.e., as de- 
scribed by Equation 7 for our data). This process creates a 
scattered galaxy number distribution given by 



S(£) = N 1 V(x)y^n(x',y,z) 



(26) 



where the normalization constant N\ can be calculated by 
requiring that S(x) = fN. Equations 25 and 26 have 
Fourier transforms 

n(k) = N[W(k) + C(k)] (27) 

S(k) = fNV(k X ) [W(0,ky,k Z ) + C(0,ky,k Z j] (28) 

where we have defined the convolved density field 

C(k) = W(k> - k) S(k') (29) 

y 

It is convenient to subtract the selection function contribu- 
tion: 



n'(fc) = n(k) - NW(k) 

S'(k) = S(k)- fNW(0,k y ,k z )V(k x ) 



(30) 
(31) 



We now construct the combined density field m(x) = (1 — 
f)n(x)+S(x) and write (for the purposes of this calculation) 
a simple power spectrum estimator 



Pcst(fc) 



lm'(fc)| 2 - N 
TP 



(32) 



where m'(k) = (1 - f)n'(k) + S'{k). The terms needed to 
determine |m'(fc)| 2 are: 



|n'(fc)| 2 = (1 - f) 2 N 2 \C(k)\ 2 

h"(k)S'(k) = f(l-f)N 2 V(k x )P'(k) 

\S'(k)\ 2 = f 2 N 2 \C(0,k y ,k z )\ 2 \V(k x )\ 2 
where 



P'{k) = Y,P{k')Wg.uWi_ 



(0,ky,kz) 



(33) 
(34) 
(35) 

(36) 



Writing the convolved power spectrum as P c (fc) = |(7(fe)| 2 , 
the final result for the power spectrum estimate is: 

Pt(fc) = (1~ f) 2 P c (k) + 2f(l- f)Re[V(k x )P'(k)] 

+ fP c (0,k y ,k z )\V(k x )\ 2 (37) 

In the special case that the selection function and scattering 
functions are constant we find that: 



Pc S t(k) = P(k) (k x = 0) 

= (l-f) 2 P(k) (k x ^0) 



(38) 
(39) 



If the Fourier modes k are binned in spherical shells as a 
function of k = |fc| then, as the value of k increases, the con- 
tribution of the mode with k x = becomes less important 
and P es t(fc) = (1 — f) 2 P(k) is an increasingly good approx- 
imation. 

We now consider the case where a fraction / of redshifts 
are scattered by a fixed displacement along the a;-axis, i.e. 
a convolution of the density field by a function U . Equation 
26 is replaced by the relation 



S(x) = JVi Y n(x',y, z)U(x - x) 



(40) 



with the Fourier transform 

S(k) = fNU(k x )[W(k) + C(k)] (41) 
Defining S'(k) = S(k) - fNU(k x )W(k) we find that 



P cst (k)/P c (k) = (l-f) 2 + 2f(l-f)Re[U(k x )} 
+ f\U(k x )\ 2 



(42) 
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For a simple shift U(x) = 8(x - 



s 



we find that 



Peat(fe) = [1 - 2/(1 - /)[1 - COS (fc^o)] P(fc) 



(43) 



These formulae are only approximations in the case of 
real data. Firstly, the sky is curved and the redshift scatters 
do not happen along a single axis of the cube. Secondly, the 
blunder fraction depends on redshift. Thirdly, the blunders 
due to line confusion are not a strict convolution of the den- 
sity field, but a transformation in redshift z\ — > Z2 of the 
form Z2 = C(l + zi) — 1, where C is a constant depending 
on the rest wavelengths of the lines. Fourthly, if a convolved 
redshift is shifted beyond the edge of the density cube, it 
does not "wrap around" as required by periodic boundary 
conditions. Fifthly, FKP estimation of the power spectrum 
is used rather than Equation 32. 

In order to measure the distortion of our measured 
power spectrum due to redshift blunders for the real data, 
we therefore created Monte Carlo simulations of galaxy cat- 
alogues with a known input power spectrum and the same 
selection function W(x) as each of our survey regions. We 
then applied the redshift blunder distributions of Figures 11 
and 12 to the mock catalogues and re-measured the power 
spectrum. Comparison of the input and output power spec- 
tra, averaging over many Monte Carlo realizations, provided 
the correction factor due to redshift blunders. We tested our 
code by reproducing the relations given in Equations 37 and 
42 in the flat-sky case. 

Figure 13 plots the correction factors for the "angle- 
averaged" power spectrum P(k) for each of the survey re- 
gions (taking a redshift interval 0.3 < z < 0.9). We see that 
a good approximation for the correction for scales k > 0.05 h 
Mpc -1 is a constant, although a more significant correction 
is required for large scales k < 0.05 h Mpc -1 . For k > 0.05 h 
Mpc -1 the approximately constant correction factor is not 
exactly equal to (1— /) -2 as predicted by Equation 39, where 
/ is the average redshift blunder rate of the catalogue; this 
is due to the mis-estimation of the denominator of Equation 
13 that occurs because W(x) is determined from the galaxy 
redshift distribution including blunders, as described in Sec- 
tion 2.5. The Monte Carlo simulations performed here also 
correct for this small bias in power spectrum estimation. 

3.3 Power spectrum measurement 

In this study we analyzed a galaxy sample drawn from Wig- 
gleZ survey observations prior to July 2009 in SDSS regions 
of our optical imaging (9-hr, 11-hr, 15-hr). Figure 14 plots 
the (R.A., Dec.) distribution of these redshifts. We imposed 
the redshift cut 0.3 < z < 0.9 in order to remove the tails of 
the redshift distribution which contain relatively few galax- 
ies. A total of TV = 56,159 galaxy redshifts remained. We 
then split the sample into three redshift slices 0.3 < z < 0.5, 
0.5 < z < 0.7 and 0.7 < z < 0.9. 

We determined the effective redshift z e g of our power 
spectrum estimate in each redshift slice by weighting each 
pixel in our 3D selection function by its contribution to the 
power spectrum error: 



z e s(k) = ^ z x 



n 9 (x)P(fc) 
l + n g (x)P(k) 



(44) 




Figure 13. The power spectrum correction factor due to redshift 
blunders for each of the survey regions analyzed in this paper, for 
a redshift range 0.3 < z < 0.9. The measured power spectrum 
must be divided by this factor in order to obtain an unbiased 
estimate of the true power spectrum. 



in each grid cell and P(k) is the power spectrum amplitude. 
In each case we used the best-fitting model power spectrum 
determined below. We evaluated this function at k — 0.15 h 
Mpc - , although the dependence on scale is weak. The ef- 
fective redshifts of each slice determined using equation 44 
are z cS = (0.42,0.59,0.78). 

We analyzed the three WiggleZ survey regions indepen- 
dently, resulting in a total of nine power spectrum mea- 
surements. We estimated the power spectrum up to a max- 
imum Fourier wavescale fc max = 0.4 h Mpc - , assuming 
the value P = 2500 ft -3 Mpc 3 for the weighting factor in 
Equation 9. This choice is motivated by our final measure- 
ment of the power spectrum amplitude presented below on 
scales k « 0.15 h Mpc - , but does not have a strong in- 
fluence on our results given that with the survey partially 
complete the measurements are limited by shot noise on 
most scales. Representative values for the other parame- 
ters in Section 3.1 are (L x ,L y ,L z ) = (600, 600, 300) ft -1 
Mpc, (n x ,n v ,n z ) = (256,256,128), V = 0.1ft -3 Gpc 3 and 
no = 5 x 10 -5 h? Mpc -3 . We combined the Fourier ampli- 
tudes in angle-averaged bins of width Ak = 0.01 h Mpc -1 . 

The nine power spectrum measurements are plotted in 
Figure 15 together with a power spectrum model derived 
using a "standard" set of cosmological parameters together 
with a prescription for redshift-space distortions. The de- 
tails of this model are described below in Section 4.1. The 
dashed and solid lines illustrate the input model, and the 
model convolved with the selection function for each region, 
respectively. The model provides an acceptable statistical fit 
to the measured power spectrum in each case. 

The corresponding nine covariance matrices dj are 
plotted in Figure 16 as a correlation coefficient 



r(i,j) = 



(45) 



where n g (x) = (N C N /V)W (x) is the galaxy number density 



Figure 16 demonstrates that the amplitude of the off- 
diagonal elements of the covariance matrices is small (note 
the choice of greyscale range). 

We also measured power spectra in wavevector bins 
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Figure 14. Greyscale map illustrating the completeness of the spectroscopic follow-up of the WiggleZ targets for the three survey regions 
analyzed in this paper. This Figure is generated by taking the ratio of the galaxy densities in the redshift and parent catalogues in small 
cells. The rr-axis and j/-axis of each panel are right ascension and declination, respectively. 



(fcperp,fcpar) perpendicular and parallel to the line-of-sight, 
respectively (we now use Fourier bins of width Afc = 0.02 h 
Mpc -1 in each direction to increase the signal-to-noise ratio 
in each bin). This 2D power spectrum allows us to recover 
the redshift-space distortion parameters which produce an 
anisotropic galaxy power spectrum. Since in our analysis we 
orient the x-axis parallel to the line-of-sight to the centre 
of each su rvey reg ion, we make the flat-sky approximation 
fcperp = y/ky + kz , fcpar = | fcjr | in this analysis. 

3.4 Systematic error study 

We investigated the dependence of our power spectrum mea- 
surement on potential systematic errors in the survey selec- 
tion function. In order to do this we re-constructed four dif- 
ferent selection functions for the 9-hr region, analyzing the 
full redshift range 0.3 < z < 0.9, with (extreme) variations 
in the method: 

• We removed the variation of the completeness of the 
parent GALEX-SDSS sample with dust and GALEX expo- 
sure time described in Section 2.2, and instead assumed a 
constant density map. 

• We removed the variation of redshift completeness 
across the field-of-view of the 2dF spectrograph described in 
Section 2.4, and instead assumed a constant redshift com- 
pleteness. 

• We removed the dependence of the survey redshift dis- 
tribution on angular position described in Section 2.5, and 
instead used a position-independent radial selection func- 
tion. 



• We parameterized the redshift distribution for each 
magnitude slice using a 9th-order polynomial rather than 
a sum of two Gaussian functions. 

In Figure 17 we plot the difference between the resulting 
power spectra measured for these four different selection 
functions and the fiducial power spectrum, in units of the 
standard deviation of the measurement at each scale. We 
note that these (extreme) variations in our understanding of 
the selection function typically cause up to « 0.5-cr shifts in 
the power spectrum estimate. The exception is the param- 
eterization of the redshift distribution, which causes large 
deviations in the very large-scale (k < 0.03 h Mpc -1 ) power 
spectrum. For the example being studied, we noted that 
the 9th-order polynomial fit was in fact providing a poorer 
match than the double Gaussian function to the observed 
redshift distributions, causing a spurious increase in mea- 
sured large-scale power. 

We conclude that our power spectrum measurements 
are likely to be robust against reasonable systematic vari- 
ations in the selection function methodology. The selection 
function in Fourier space is compact and the corrections are 
only significant on the largest scales. 



4 COSMOLOGICAL PARAMETER FITS 

In this Section we present some initial comparisons between 
our power spectrum measurements and cosmological models. 
Our aim is to establish whether or not our measurements 
are consistent with the predictions of the standard ACDM 
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Figure 15. The ID galaxy power spectrum for each of the survey regions and redshift slices analyzed in this paper. The data points in 
each plot are the power spectrum measured by FKP estimation. The dashed line is an (unconvolved) model power spectrum evaluated 
by integrating equation 46 over angles. We use cosmological parameters f2 m = 0.3, f2t,/f2 m = 0.166, h = 0.72, n B = 0.96 and erg = 0.8, 
together with values of the redshift-space distortion parameters and linear bias factor fit to the 2D power spectrum for each redshift 
slice. The solid line is the result of convolving this model power spectrum with the selection function for each region. The chi-squared 
statistic is calculated over the range k < 0.4 h Mpc -1 using the full covariance matrix. The x-axis and y-axis of each panel are Fourier 
wavescale k in units of h Mpc -1 and power spectrum amplitude P(k) in units of h~ 3 Mpc 3 , respectively. 
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Figure 16. Greyscale plot of the correlation coefficient r of Equation 45 for each of the survey regions and redshift slices analyzed in 
this paper, indicating the degree of covariance between the power spectrum measurement in different Fourier bins. Note the choice of 
greyscale range such that the darkest shade corresponds to r = 0.2. 
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Figure 17. Variations in the power spectrum estimate for the 
9-hr region, in units of standard deviations in the fiducial power 
spectrum, obtained by changes in the methodology used to con- 
struct the selection function. 

framework with parameters determined by observations of 
the Cosmic Microwave Background radiation. Future papers 
will perform a more comprehensive analysis. 

4.1 Model power spectra 

We derived model matter power spectra Pm(k) as a func- 
tion of redshift using the CAMB software package (Lewis, 
Challinor & Lasenby 2000) which is based on CMBFAST (Sel- 
jak & Zaldarriaga 1996), including corrections for non-linear 
growth of structure using the fitting formulae of Smith et al. 
(2003) (halofit=l in CAMB). Throughout this analysis (un- 
less otherwise noted) we fix values for the Hubble parameter 
h = 0.72, scalar index of spectral fluctuations n B = 0.96 and 
normalization <rg = 0.8, and vary only the matter density fi m 
and the baryon fraction flb/fi m . Our choices of parameter 
values are consistent with recent observations of the Cosmic 
Microwave Background radiation (Komatsu et al. 2009). 

For this initial parameter fitting we assumed a linear 
bias factor b for the galaxy population. The assumption 
of linear bias will be investigated in more detail in future 
studies through the use of mock catalogues generated from 
dark matter simulations, analysis of the galaxy bispectrum 
(e.g. Verde et al. 2002; Nishimichi et al. 2007), and mea- 
surements of the halo occupation distribution of WiggleZ 
galaxies (e.g. Zehavi et al. 2005; Blake, Collister & Lahav 
2008) . We note that in the current analysis the value of the 
bias factor is significantly degenerate with the choice of as 
(such that b oc l/as). 

We modified this real-space matter power spectrum to 
include redshift-space distortions, which we characterized in 
the standard manner using two parameters. The first pa- 
rameter, /, models the effect of large-scale coherent infall 
velocities. In linear theory for the growth of fluctuations, / 
is related to the linear growth factor D(a) at scale factor a by 
/ = din D/d\na. This functional form is well-approximated 
by / = S7 m (z) 7 where Q m (z) is the matter density measured 
by an observer at redshift z, and 7 as 0.55 for a standard 
ACDM cosmology. Coherent velocities can also be param- 
eterized by /3 — f/b, introducing a dependence on galaxy 



type via the bias factor. The second parameter, u v (in units 
of h km s _1 ), describes small-scale random virialized veloci- 
ties, for which we assumed an exponential pairwise distribu- 
tion. The overall effect on the galaxy power spectrum P ga i is 
to induce an anisotropic distortion dependent on the angle 
9 of the wavevector to the line-of-sight, parameterized by 
H — cos 6: 

P ^ = b2p ^ i%Hoi%)' (46) 

4.2 Fits to the 2D power spectra in redshift slices 

We fitted these models to the nine independent 2D power 
spectra split into tangential and radial components, in order 
to determine the redshift-space distortion parameters and 
test the ACDM self-consistency relation /3 « f2 m (z) - 55 /b 
as a function of redshift. In this initial analysis we fixed 
the matter density parameter f2 m = 0.3 and baryon frac- 
tion S7b/f2 m = 0.166 (similar to the parameters listed in 
Komatsu et al. (2009) that provide the best fit to the 5- 
year data set of the Wilkinson Microwave Anisotropy Probe 
(WMAP); we fix f2 m = 0.3 to match the value assumed for 
the distance-redshift relation in our power spectrum anal- 
ysis). We then varied the redshift-space parameters / and 
g v and the bias factor b, convolving each model power spec- 
trum with the survey selection function for each region and 
minimizing the chi-squared statistic calculated using the full 
covariance matrices. For each redshift slice we combined the 
results from the three survey regions. The best-fitting mod- 
els produce a good fit to the measured power spectrum, as 
evidenced by the values of the chi-squared statistic in the 
three redshift slices: X 2 = (1113.9,1229.6,1069.2) for 1137 
degrees of freedom. 

The measurements of the growth rates at each effective 
redshift z cff = (0.42, 0.59, 0.78) are / = (0.73 ± 0.09, 0.75 ± 
0.09,0.71 ± 0.14), marginalizing over the parameters a v 
and b, where 1-a error ranges are indicated. These should 
be compared to the ACDM prediction for S7 m (0) = 0.3: 
/ = (0.72,0.78,0.83). We find that the model successfully 
passes this test of self-consistency. Our measurements of 
the growth rate are compared with the model in Figure 
18 by plotting the quantity / x as, mass (z) versus redshift. 
Previously-published results are also indicated, as summa- 
rized by Song & Percival (2009). The WiggleZ Survey data 
provides a growth measurement across the redshift range 
0.4 < z < 0.8 with a precision similar to previous work at 
z < 0.4. We highlight in particular the 20% accuracy of our 
measurement at z = 0.78. A future study will perform a 
more detailed analysis of these results. 

The fitted values of the pairwise velocity dispersion in 
each redshift slice are a v = (354 ± 41, 294 ± 31, 216 ± 58) h 
km s _1 (marginalizing over / and b). We find evidence 
that the pairwise small-scale velocity dispersion of Wig- 
gleZ galaxies systematically decreases with increasing red- 
shift. The measured bias factors for each redshift slice are 
b = (0.93 ± 0.03, 1.08 ± 0.03, 1.20 ± 0.06) (marginalizing 
over / and a v ). We can compare these galaxy bias mea- 
surements to those deduced from the WiggleZ survey small- 
scale correlation function measurements: for three redshift 
slices similar to those analyzed here, Blake et al. (2009) ob- 
tained b — (1.01, 1.27, 1.27) (although assuming a higher 
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Figure 18. Measurements of the growth rate of structure, /, ob- 
tained in three redshift slices by fitting WiggleZ survey data in 
three survey regions. The values of / are weighted by o-g mass (z) 
as discussed by Song & Percival (2009). The measurements are 
compared to results previously published for the 2dFGRS, SDSS 
LRG, and VIRMOS-VLT Deep Survey (VVDS) samples, as col- 
lected by Song & Percival. A prediction for the ACDM cosmolog- 
ical model with f2 m = 0.3 and erg = 0.8 is overplotted. 

value of erg = 0.9). The main cause of this difference is the 
superior determination of ft in the current analysis. Other 
issues wth this comparison include: (i) galaxy bias is a scale- 
dependent function, (ii) small-scale pairwise velocities were 
not been modelled in the Blake et al. (2009) analysis, (iii) 
the power-law correlation function model assumed in Blake 
et al. (2009) breaks down at large scales. 

4.3 Stacked ID power spectra 

Using these model fits we can combine the nine ID power 
spectrum measurements in the different regions and redshift 
slices into a single "stacked" survey ID power spectrum us- 
ing inverse-variance weighting. The difficulty with this step 
is that the power spectrum amplitude in each region is mod- 
ulated by a different level of convolution with the selection 
function, as illustrated by Figure 15. In addition, the shape 
of the power spectrum varies with redshift in accordance 
with the differing redshift-space distortion parameters and 
galaxy bias factors. 

Therefore, before combining the results in the different 
regions and redshift slices, we performed an approximate 
"de-convolution" of amplitude by multiplying the power 
spectra by a correction factor equal to the scale-dependent 
ratio of the convolved and unconvolved model power spectra, 
i.e. the ratio of the solid and dashed curves plotted in Fig- 
ure 15. (In detail we did not use the model power spectra to 
generate these corrections but a "reference" power spectrum 
in which the baryon oscillation features have been smoothed 
out, i.e. the "no- wiggles" power spectrum of Eisenstein & Hu 
(1998). We preferred to correct the amplitude using a wiggle- 
free power spectrum to avoid spuriously introducing appar- 
ent baryon oscillations into the combined power spectrum). 
We additionally corrected the measurement in each Fourier 
bin by the angle-averaged ratio of the redshift-space galaxy 
power spectrum (using the best-fitting values of /3 and a v for 



each redshift slice) to the real-space matter power spectrum 
at redshift z — 0.6. Following these corrections, the nine 
power spectra "line up" with consistent shape and ampli- 
tude, and can be combined. In Figure 19 we show separate 
combinations of power spectra across survey regions in each 
redshift slice, and across redshift slices in each survey re- 
gion. The combination of all nine power spectra is presented 
in Figure 20; the lower panel indicates the fractional accu- 
racy of the measurement, which approaches 5% in Fourier 
bins of width Afc = 0.01 h Mpc - . The scale dependence of 
the fractional accuracy is determined by a balance between 
the increasing number of Fourier modes contributing to each 
successive bin (tending to decrease the error), and the in- 
creasing importance with k of shot noise relative to cosmic 
variance (tending to increase the error). Finally, in Figure 
21 we display the 2D power spectra for each redshift slice, 
combining different regions in the manner described above. 

4.4 Fits for matter and baryon densities 

As a final exercise we fitted the measured power spectra 
for the cosmological matter and baryon densities (parame- 
terized by fi m and /b = Qb/Q m ). We tried three different 
initial approaches to this analysis, with increasing degrees 
of sophistication: 

• First, we generated unconvolved, real-space, non-linear 
power spectra and fitted the combined survey data points 
plotted in the upper panel of Figure 20 by varying fl m , /b 
and the bias factor b. 

• Second, we fitted the convolved, redshift-space ID 
power spectra in each of the survey regions and redshift 
slices, fixing the redshift-space distortion parameters /3 and 
a v at the best-fitting measurements using the 2D power 
spectrum fits at each redshift and varying fl m , /b and b. 

• Third, we fitted the convolved, redshift-space 2D power 
spectra, fixing the value of a v at the best-fitting values. We 
varied fl m , /b and the bias factor 6, determining the value 
P = Q m (z)°- 56 /b in each case (assuming the ACDM model). 
For each different value of f2 m , we scaled <rg assuming the 
"WMAP normalization" (Komatsu 2009, Section 5.5) and 
also scaled the geometry of the survey box in accordance 
with each trial cosmology: given a fiducial co-moving dis- 
tance Dad and Hubble parameter H^d for a redshift slice 
calculated using fi m = 0.3, and also given trial values D, 
H corresponding to a different value of fi m , the width of 
Fourier bins tangential and radial to the line-of-sight scales 
as D/Dad and Had/H, respectively; furthermore the vol- 
ume of the box (hence amplitude of the measured power 
spectrum) scales as (D/Dad) 2 x (Had/H). 

Figure 22 displays the probability contours in the 2D 
space of fi m and /b = fib/^m, marginalizing over the bias 
factor, for each of the three approaches described above. The 
different methods produce broadly consistent results, and 
the best-fitting values are in good agreement with those de- 
termined from the latest measurements of the temperature 
fluctuations of the Cosmic Microwave Background radiation 
(Komatsu et al. 2009). 

The probability contours in Figure 22 display evidence 
of the well-known degeneracy between Sl m and flb/Qm in 
the determination of the overall shape of the matter power 
spectrum, which is exacerbated by the fact that we cannot 
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Figure 19. The power spectra determined by combining the measurements in the different survey regions (for each redshift slice) and 
different redshift slices (for each survey region). Before combination, we adjust the shapes of the power spectra to allow for the differing 
degrees of convolution with the selection function and the differing redshift-space effects. The solid line in the panel has been generated 
by combining the best-fitting model power spectra in a similar manner. The x-axis and y-axis of each panel are Fourier wavescale k in 
units of h Mpc -1 and power spectrum amplitude P(k) in units of h~ 3 Mpc 3 , respectively. 
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Figure 20. Upper panel: The power spectrum determined by combining the measurements in the different survey regions and redshift 
slices analyzed in this paper using inverse-variance weighting. Before combination, we adjust the shapes of the power spectra to allow 
for the differing degrees of convolution with the selection function and the differing redshift-space effects. The solid line in the panel has 
been generated by combining the best-fitting model power spectra in a similar manner. Lower panel: The fractional error in the combined 
power spectrum as a function of Fourier wavenumber k. 
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Figure 21. The 2D galaxy power spectra for each of the redshift slices analyzed in this paper, obtained by combining measurements 
in different survey regions, and plotted as a function of wavevectors (k peTp , fc par ) tangential and radial to the line-of-sight, respectively. 
The function is represented using both greyscale and contours; the contour levels are logarithmically spaced between P = 1000 and 
10,000ft- 3 Mpc 3 . The non-circularity of the contours encodes the imprint of large-scale galaxy peculiar velocities, as discussed in the 
text. 
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Figure 22. Probability contours of f2 m and Q^,/Cl m fitting 
to WiggleZ survey power spectra using the three different ap- 
proaches described in the text. The inner and outer contours for 
each set enclose 68% and 95% of the likelihood, respectively. 



yet detect the imprint of the baryon acoustic oscillations in 
our power spectrum measurement. This degeneracy will be 
broken as the WiggleZ Survey progresses. 



5 CONCLUSIONS 

In this paper we have described our method of determining 
the selection function of the WiggleZ Dark Energy Survey, 
and have presented the current measurement of the large- 
scale galaxy power spectrum using 56,159 redshifts of bright 
emission-line galaxies spanning redshifts 0.3 < z < 0.9. 
This sample constitutes approximately 25% of the final Wig- 
gleZ survey. We have quantified and categorized the redshift 
blunder rate and determined its effect on the power spec- 
trum measurement via analytical calculations and detailed 
simulations. We conclude that: 

• The selection function of the WiggleZ survey is compli- 
cated by the proximity of the faint magnitude threshold to 
the completeness limit of the input catalogues, in particular 
for the GALEX UV data. We quantified the incompleteness 
in the parent target catalogue as a function of GALEX ex- 
posure time and Galactic extinction via fitting formulae. 

• We adopted a Monte Carlo technique to determine the 
relative completeness of the spectroscopic follow-up at any 
position. This technique allows for the complex overlapping 
of survey pointings and for the systematic variation of red- 
shift completeness across the 2-degree field-of-view of the in- 
strument. We also allowed for the magnitude prioritization 
of the spectroscopic follow-up which results in a position- 
dependent galaxy redshift distribution. 

• The WiggleZ survey contains redshift blunders resulting 
from emission-line confusion (most significantly, [OIII], H/3 
and Ha mis- identified as [Oil]) and from sky emission lines 
mis-identified as [Oil]. The overall blunder rate is about 5%. 
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The effect of the redshift blunders on the power spectrum 
measurement is well-approximated as a constant reduction 
in amplitude for scales k > 0.05 h Mpc" 1 combined with 
an enhanced level of reduction for large scales k < 0.05 h 
Mpc" 1 . 

• We measured ID (angle-averaged) and 2D (binned in 
tangential and radial modes) galaxy power spectra for nine 
independent survey regions and redshift slices using the 
method of Feldman, Kaiser & Peacock (1994). The ID power 
spectrum for the whole sample, combining these measure- 
ments, has a fractional accuracy of about 5% in Fourier 
bins of width Ak = 0.01 h Mpc -1 . The 2D power spectra 
show the expected anisotropic signatures of redshift-space 
distortions due to large-scale coherent infall and small-scale 
virialized motions. 

• The power spectrum data are well-described by a model 
power spectrum with matter and baryon densities consis- 
tent with those determined from observations of the Cos- 
mic Microwave Background radiation. The model includes 
non-linear corrections, redshift-space distortions and a lin- 
ear galaxy bias factor. 

• The 2D power spectra allow us to measure the growth 
rate of structure across the redshift range 0.4 < z < 0.8. We 
obtain results similar in precision to previous determinations 
at z < 0.4, including a measurement at z = 0.78 with 20% 
accuracy. 

Future studies will present full cosmological parameter fits 
and combinations of these results with other datasets, in- 
cluding implications for the growth of cosmic structure and 
Gaussianity of the initial conditions, and extend these anal- 
yses to the final WiggleZ survey catalogues. 
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